Phenomenology of the Equivalence Principle with Light Scalars 
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Light scalar particles with couplings of sub-gravitational strength, which can generically be called 
'dilatons', can produce violations of the equivalence principle. However, in order to understand 
experimental sensitivities one must know the coupling of these scalars to atomic systems. We report 
here on a study of the required couplings. We give a general Lagrangian with five independent 
dilaton parameters and calculate the "dilaton charge" of atomic systems for each of these. Two 
combinations are particularly important. One is due to the variations in the nuclear binding energy, 
with a sensitivity scaling with the atomic number as A^ 1 ' 3 . The other is due to electromagnetism. 
We compare limits on the dilaton parameters from existing experiments. 



The Equivalence Principle (EP) is one of the most 
exquisitely tested features in science, with a present sen- 
sitivity of 
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for comparisons of the acceleration of two test masses in 
a gravitational field [l|, |2(. Further tests of this prin- 
ciple remain important because of this high precision 
- the tests can be sensitive to very small effects that 
other experiments might not notice [3, |j| . Indeed, new 
space-based experiments such as MICROSCOPE Q, the 
Galileo Galilei project |6| and STEP Ml as well as inno- 
vative experiments with cold atoms [a [9( and sub-orbital 
rockets [10( are being planned that could push the sensi- 
tivity several orders of magnitude further. 

In comparing experiments, it is common to invoke a 
hypothetical model with the exchange of a light vector 
particle coupled to baryon number. This is convenient 
because it is easy to know the baryon number of a given 
material, and hence the coupling of this vector particle 
to atomic systems is known. However this is an unlikely 
physical situation, as very light and very weakly cou- 
pled vectors are rare in modern theories. A more plau- 
sible candidate as a source of EP violation is a light 1 
scalar field with a coupling to matter that is weaker than 
gravitational strength. We will refer to these generically 
as 'dilatons', although they may have origins other than 
string theory or models involving dilation symmetry [lll - 
[l4| . Being Lorentz scalars, dilatons will couple to Lorentz 
invariant combinations of the Standard Model fields, the 
quarks, gluons and photons. Significant theoretical effort 
is needed to translate these fundamental interactions into 
the required couplings of the atomic systems used in ex- 



1 We will generally assume in the following that the scalar field 
we consider is essentially massless on the scales that we discuss, 
although is it simple to modify this assumption. 



periments. In [15[ we provide an extensive study of the 
theoretical ingredients to this program, and in this pa- 
per we summarize the results which are most useful to 
experiments. 

The important physics is readily understood, even if it 
is somewhat more difficult to calculate. Most of the mass 
of atomic systems comes from the gluonic QCD interac- 
tions, which sets the scale for the strong interactions. 
However, atomic masses also have smaller dependencies 
on the masses of the light quarks and on the electric 
charge. This can occur either in the masses of the indi- 
vidual nucleons or in the binding energy of the nucleus. 
General relativity couples universally to all forms of en- 
ergy, but dilaton couplings can and will vary among the 
different ingredients to the total atomic mass. We find 
particular sensitivity in the binding energy coming from 
the quark mass dependence, and this is the most novel 
aspect of our analysis 2 . As described below, in addi- 
tion to the well known dependence on electromagnetism, 
the binding energy has an enhanced dependence on the 
quark masses because of the importance of the pion in 
nuclear binding and the strong dependence of the pion 
mass on the light quark masses. Although we provide a 
description of all effects, we will also provide a simpler 
parameterization of the dominant hadronic binding sen- 
sitivity, proportional to A -1 / 3 , and the electromagnetic 
sensitivy, proportional to Z(Z — 1)/A 4 / 3 . 

A massless dilaton <j> modifies the Newtonian interac- 
tion between a mass .4 and a mass B, into the form (see. 
e.g. ET 



V = -& (1 + a A a B )- 

TAB 



(2) 



If the dilaton mass is important the second term includes 
an extra exponential factor exp(— th^tab)- In this inter- 



2 Damour Q| and Dent Il6( have highlighted this need for the study 
of the nuclear binding energies. 



action potential, the scalar coupling to matter is mea- 
sured by the dimensionless factors, a^g. To calculate 
these, we add the dilaton Lagrangian to that of the Stan- 
dard Model and use the result to evaluate the total mass 
of the atomic system. The a a coupling is then found via 
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Here, K = \/4nG is the inverse of the Planck mass 3 so 
that the products ktua and kc(> are dimensionless. This 
ensures that this definition of a a is valid in any choice of 
units. In terms of the ct^'s, the violation of the (weak) 
EP, i.e. the fractional difference between the accelera- 
tions of two bodies A and B falling in the gravitational 
field generated by an external body E, reads 
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In the last (approximate) equation we have assumed that 
the a's are small. 

We couple the dilaton to the light fields of the Stan- 
dard Model. The heavy quarks and weak gauge bosons 
are assumed to be integrated out and we assume that 
the dilaton theory has been matched to the light fields 
below the scale of the heavy quarks. This procedure in- 
troduces five dimensionless dilaton-coupling parameters, 
d e ,d g for the couplings to the electromagnetic and glu- 
onic field-strength terms, and d m< , , d mu , d md for the cou- 
plings to the fermionic mass terms. [We are using here 
the fact that a (^—dependent coupling to the kinetic term 
of a fermion, f(cf>)ipilfiip, can be absorbed in a suitable 
</>— dependent rescaling of tp.] We add these interac- 
tions to that of the Standard Model, C — Csm + Ant 
and normalize these five dimensionless dilaton parame- 
ters d e ,dg,d mc ,d mu ,d md so that they correspond (when 
considering the linear couplings to (f>) to 
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These are chosen to correspond to renormalization group 
invariants [lj, [l5|. For the quark mass terms, physics 
tells us that it is preferable to work with the symmetric 
and antisymmetric combinations 



m = -(m d + m u ) 



5m = (m d - m u ) (6) 



We use units such that 



so that we define dilaton parameters dm, d$ m conjugate 
to these combinations 



Ch 



d m rh(dd + uu) H — —5m(dd — uu) 



(7) 
and use these parameters instead of d mu , d md . 

We need to calculate the effect of these couplings for 
atoms. In order to do this we transform the scalar field 
dependence into an implicit dependence on the param- 
eters of the Standard Model, and use the decades of 
research connecting the Standard Model to observable 
physics. For example, we follow the pioneering method 
of [ll| for the electromagnetic coupling 
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where the last equality is valid at the linear level in n<j) 
(which is the level at which we define the dilaton cou- 
plings here) . As we work with a rescaled electromagnetic 
field (A hcic = eA usual ), the only location where the elec- 
tric charge occurs in the Lagrangian is the one explicitly 
shown above. This allows the dilaton field to be absorbed 
into the following <f> dependence of the fine-structure con- 
stant a — e 2 /(47r) 



a(4>) 



(1 + d e K(b)a 



(9) 



and allows us to obtain the dilaton coupling from the 
dependence of the masses on the electric charge 



„(d e 



d e dm,A 
rriA da 



(10) 



Likewise the fermion mass terms are normalized such 
that we can define a field-dependent mass 



mi(4>) = (1 + d mi K(j))m l 



u.d.e 



(11) 



and we obtain the dilaton coupling from the dependence 
of the atomic masses on the fermion masses. Finally, the 
terms proportional to the gluonic coupling d g are normal- 
ized so that they are proportional to the QCD trace of 
the energy momentum tensor [1J, [l5[ , which allows the 
matrix element to be readily calculated. This turns the 
parameter d g into one that measures the <f> sensitivity of 
the QCD mass scale, say A3. We then have 



9 In A 3 
d[K</>] 
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(12) 



Variation of dimensionless ratios such as ln(mi/A3) then 
involve differences of the coupling parameters d m — d g . 

The structure of the coupling to individual nucleons 
is reasonably well known. The quark mass contributions 
are known from the nucleon sigma term [17[ and from 
the baryon mass splittings, and electromagnetic contri- 
butions to the masses have been estimated [18l ]. 



Nuclear binding is well parameterized by the semi- 
empirical mass formula 4 , 
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Typical fit values for these parameters are [19j a v = 
16 MeV, a s = 17 MeV, a a = 23 MeV, a p = 12 MeV, a c = 
0.717 MeV. The Coulomb term a c is linear in the electro- 
magnetic fine structure constant and therefore directly 
yields the electromagnetic coupling that we need. For 
the hadronic component, we use our previous work on 
the quark-mass dependence of nuclear binding [20j . The 
primary physics here is that the central nuclear potential 
has an important component from the exchange of two 
pions. Because of the nature of chiral symmetry in QCD, 
the pion mass-squared is directly proportional to the av- 
erage light quark mass, m^ ~ rh. In the overall binding 
energy, there is a partial cancelation between the attrac- 
tive two pion component, which is particularly sensitive 
to the pion mass, and the repulsive short range potential, 
which is less sensitive. The variation with pion mass then 
is stronger than simple expectations. We have estimated 
the variation of both the central terms a v , a s and the 
asymmetry energy a a and find the following contribution 
to Ola 
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(14) 



Here we use the notation 



m A 



(15) 



where we take m amu = 931 MeV as the nucleon mass with 
the average binding energy, 8 MeV, subtracted. This 
factor is very close to unity throughout the periodic table. 
We estimate a 30% uncertainty in these numbers, and 
note that our formalism in [15| is most reliable for heavier 
elements. 

In order to highlight the difference between various 
materials we define 



a A = d q 



OLA 



The overall common coupling d g does not violate the EP, 
while the EP variation is contained in a. The result of our 
calculation can be summarized in four 'dilaton charges' 
Qm, Qsm, Qm c , Qe- Each charge gives the strength 
of the EP-violating coupling corresponding to a given 
dilaton parameter for an atom of charge Z and atomic 



4 For simplicity in this Letter we drop the pairing interaction as it 
does not play a significant role in our analysis 



number A, namely 

OLA = [{dm - d g )Qm + {dsm - d g )Qsm 

The four dilaton charges are given by 



(17) 



= F A 
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and 



Qe = F A 



-1.4 + 8.21 + 7.7^1 



x 10" 



(18) 



(19) 



(20) 



(21) 



Again note that Fa (see Eq. (|15p ) can be readily approx- 
imated as unity. 

Inspection of our general results of the previous para- 
graph reveals that many terms have small effects. In- 
deed, insertion of representative values of A, Z indicate 
that there are two dominant effects - those of the A^ 1 ^ 3 
and Z(Z — 1)/A 4 ' 3 variation in the binding energy, from 
strong and electromagnetic effects respectively. Until 
such fortunate time when we need high precision to ac- 
count for multiple measurements, we can obtain a more 
useful formula by truncating to these two terms. Doing 
so yields a simpler, approximate formula involving only 
two dilaton charges 



OLA 



d* + [(d m -d g )Q[ h + d e Q' e ] 



elA 



where 



d* g =d g + 0.093(d A - d g ) + 0.000274 

and where 
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(16) and 
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(22) 
(23) 

(24) 
(25) 



We think that these approximate expressions capture all 
the potentially dominant EP violation effects. These 
dilaton charges for many materials are shown in Table 
I. 

The signals for EP-violation then has two terms, 
namely 

— J = (olb - OL C )oL E = [D m Q' A + D e Q' e ] BC (26) 

a / BC 



TABLE I: Approximate EP-violating 'dilaton charges' for a 
sample of materials. These charges are averaged over the 
(isotopic or chemical, for SiCh) composition. 
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^Crh 


Q'e 


Li 


7 


3 


18.88 xlO" 3 


0.345 xlO" 3 


Be 


9 


4 


17.40 xlCT 3 


0.494 xlO" 3 


Al 


27 


13 


12.27 xlO" 3 


1.48 xlO" 3 


Si 


28.1 


14 


12.1 xlO" 3 


1.64 xlO" 3 


Si0 2 






13.39 xlO" 3 


1.34 xlO" 3 


Ti 


47.9 


22 


10.28 xlO" 3 


2.04 xlO" 3 


Fe 


56 


26 


9.83 xlO" 3 


2.34 xlO" 3 


Cu 


63.6 


29 


9.47 xlO" 3 


2.46 xlO" 3 


Cs 


133 


55 


7.67 xlO" 3 


3.37 xlO" 3 


Ft 


195.1 


78 


6.95 xlO" 3 


4.09 xlO" 3 



where [Q]sc = Qb — Qc, and where the ' dilaton charges' 
are (approximately) given by Eq. (|24|) and Eg. ((25)) . The 
corresponding dilaton coefficients Di are given by 



Uir, 



= d* (dra - d g ), 



D, 



= d* g d e 



(27) 



If we were assuming that the dilaton parameter d e is 
much smaller than dm — d g , we could go further and 
conclude (in view of the numerical results indicated in 
Table I) that the signal Q' e is sub-dominant w.r.t. Q m . 
In that case we would end up with a uni-dimcnsional EP 
signal proportional to [Q' m ]BC- 

As an application of our results, let us extract the con- 
straints on the dilaton parameters from the most sensi- 
tive present experiments. The EotWash experiment |l| 
compares Be (A=9, Z=4) and Ti (A=47.9, Z=22), with 
the constraint 

— ) = (u Be - a Tl )a Eai th = (0.3 ± 1.8) X 10" 13 

a /BcTi 

(28) 

Working at the two-sigma level, i.e. (0.3 ± 3.6) x 10 13 , 
and neglecting the central value 0.3, the rewriting of this 
equation in terms of the theoretical dilaton parameters 
D m ,D e yields 



\D m + 0.22D e \ < 5.1 x 10 



-it 



(29) 



The Lunar Laser Ranging experiment [2J compares the 
acceleration of the Earth and the Moon towards the Sun, 



with the constraint 

= (^Earth — aMoon)aSu 
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(-1.0±1.4)xl0 



-13 



LLR 



(30) 

We approximate the Earth's mantle composition as being 
SiC>2, and the Earth'core as being iron, with the Moon 
being similar to the mantle. The constraint here is 



\D m + 0.28D e \ < 24.6 x 10 



-li 



(31) 



We see that although the Lunar experiment has a 
slightly better differential-acceleration sensitivity, the 
laboratory-based test is more sensitive to the dilaton co- 
efficients because of a greater difference in the dilaton 
charges of the materials used, and of the fact that only 
one-third of the Earth mass is made of a different mate- 
rial. Note that, in the dilaton models considered here, the 
EP-violation associated to the gravitational self energy 



21 



is negligible compared to the the matter couplings 
that we have calculated above. 

In summary we have given a general description of EP 
violations associated to the possible couplings of light 
scalars, and in particular have calculated, for the first 
time, the effects of nuclear binding (using recent work 
on its quark-mass dependence). The resulting A" 1 / 3 de- 
pendence is seen to be one of the major factors for EP 
violations due to the dilaton couplings, along with the 
previously quantified Z(Z — l)/^4 4 / 3 dependence of the 
electromagnetic interaction. We expect that our results 
will be useful in assessing the optimal choices for ma- 
terials in future experiments, and in interpreting their 
results. These upcoming experiments have the poten- 
tial to probe new territory in fundamental physics and 
may have the opportunity to uncover rare and exciting 
physics. 
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